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We study the strong correlation effects in the vicinity of the Mott metal- insulator transition using 
^rr ^ ' coupled clean or disordered Hubbard chains with a infinitely large coordinate number D±_ — > oo in 

the direction perpendicular to the chains and with a long-range transverse hopping. Strong electron 
correlation effects are treated partially non-perturbatively with the use of the exact results for the 
ID Hubbard model. In the case of clean systems, the thermodynamic and transport quantities which 
characterize the Mott transition from the Fermi liquid state, such as the specific heat coefficient, 
the Drude weight, and the compressibility, are obtained as functions of hole-doping 5 = 1 — n (n, 
electron density) by the systematic expansion in terms of the inverse of the transverse hopping range 
l±. We find that the (5-dependence of these quantities shows non- universal behaviors with exponents 
depending on the strength of electron-electron interaction. In the presence of disorder, it is shown 
that the frequency-dependence of the dynamical conductivity obeys the Mott's law, implying the 
flj , possibility of Mott glass state for 5^0, provided that there exists a finite range interaction. 
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I. INTRODUCTION 



o 

^^ in general, the concrete computation of these quantities in the vicinity of the MH transition is not tractable. It 
^z" difficult to go beyond mean field theories such as the Gutzwiller-type approximation in dimensions larger than one 



The Mott- Hubbard (MH) metal- insulator transition is a long-standing important issue in strongly correlated electron 
systems JJq There, has been much progress in the theoretical studjf of the MH transition in the one-dimcnsional(lD) 
' O Hubbard modelBia and the infinite-dimensional Hubbard model. hM However, our understanding for the MH transition 
5 in the spatial dimensions between one and infinity is still far from complete. In the metallic state near the MH 

transition, it is expected that the tendency toward the MH transition may appear in thermodynamic and transport 
properties such as specific heat coefficient, compressibility, and Drude weight in the metallic state, as the electron 
density approaches to the half-filling, or the strength of electron correlation approaches to a critical value Uc- However, 
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The purpose of this paper is to propose and solve the model in which an advanced treatment of electron correlation 
^5 ] beyond conventional mean field theories is possible in dimensions larger than one. 
f^ , Here, we consider the coupled Hubbard chains with a large coordinate number 2D±^ in the direction perpendicular 

T-H ■ to the chains, and with a long-range transverse hopping. We investigate how the sign of the MH transition manifests 
0^ , in highly anisotropic Fermi liquids, as the electron density n approaches the half- filling. Our strategy is as follows. 
^^ ■ We expand the partition function in terms of the transverse hopping. Then, we can deal with the intra-chain 

electron-electron scattering processes exactly exploiting the Bethe ansatz exact solution for the ID Hubbard model. 

The inter-chain electron-electron scattering processes are treated perturbatively. This partially exact treatment of 
G ' electron-electron interaction enable us to study strong correlation effects in the vicinity of the MH transition beyond 
1-^ , simple mean field theories, 
rj ' We consider the following model Hamiltonian, 

O '• 

O I H = -t^4ijCa»+ij + h.c. + uY^c\^jC'^i^jc\i^^cii,j, 

[X^ i -Ylt^JJ'^'UjCai,!'- (1) 

Here we introduce the long-range transverse hopping with a hopping range l± , 

^^(^-l) = 1 — m= XI y^cosnk^j. (2) 

n=-l± j = l 

Several years ago, Boies et al. studied the two-dimensional version of the above model in the context of dimensional 
crossover, and pointed out that in the limit of infinitely long-range transverse hopping l± —^ oo, the mean field 
solution of the Fermi liquidi state is exact, and that the system is reduced to the non-interacting fermion system 
with a renormalized massJ13 Here we consider the expansion in l/l± around this non-interacting limit. As will be 



shown in this paper, in order to make the expansion in terms of l/l±_ controllable, we need to introduce a sufRciently 
large coordinate number in the transverse direction 2D^^. The critical dimension of D±^ will be diseussed later. 
For convenience, we consider the infinitely large D±^ which simplifies perturbative calculations in 1/^j^.ll^ The l/l±r 
corrections generate multi-particle interactions, which give rise various interesting electron correlation effects. Here, 
we consider the paramagnetic Mott transition, since any magnetic phase transitions may conceal generic effects of 
electron correlation which drive the system to the Mott insulator. 

Another purpose of this paper is to investigate effects of disorder in the vicinity of the Mott transition. Very 
recently, it is proposed by Orignac et al. that in the presence of both disorder and umklapp scattering due to electron 
correlation, a novel state called Mott glass may realizcllJ In this state, a single-particle excitation has a charge gap, 
whereas particle-hole excitation is gapless, and the dynamical conductivity shows the Mott's law, cr(a') ~ w^ for small 
to. Orignac et al. used a ID interacting spinless fermion model with disorder to discuss this possibility. However, as 
is well-knawn, the excitation gap of the ID spinless fermion model at half-filling is destroyed by infinitesimally small 
disorder.t3 To avoid this difficulty, they add the mass term to the model to stabilize the gapped state. In contrast 
to the spinlpiS fermion model, the Mott-Hubbard gap of the ID Hubbard model is stable against sufficiently small 
disorder .E3T3 Thus introducing random potential to our model, eq.(^, we can investigate more directly whether the 
Mott glass is realized or not by disorder. We would like to address this issue in this paper. 

The organization of this paper is as follows. In §2, we introduce the effective action which is the basis of our 
approach, and discuss the case of the infinitely long range transverse hopping, where the system is reduced to the 
non-interacting fermions with an enhanced mass. In §3, we consider 1/^j^-corrections to the results obtained in §2. 
We see that the multi-particle interactions generated by the I/Z^-corrections change the doping-dependence of the 
thermodynamic and transport properties which characterize the MH transition. In §4, we discuss effects of disorder 
in the vicinity of the Mott transition, and a possibility of the Mott glass state. We give a summary in §5. 

II. EFFECTIVE ACTION AND THE FERMI LIQUID STATE IN THE LIMIT OF INFINITELY 

LONG-RANGE TRANSVERSE HOPPING 

A. Effective action 

In this subsection, we give the effective action of the model (|l|) which is suitable for the following argument. In order 
to exploit the exact results of the ID Hubbard model, we expand the action in terms of the transverse hopping. This 
approach was taken before in the study of dimensional crossover between the Luttinger liquid and the Fermi liquidOa 
Here we just orrtiine the basic idea. We introduce the Grassmann variable ipaij , applying the Hubbard-Stratonovich 
transformation,^ 

Then averaging over Caij, c^j with respect to the action for the ID Hubbard model, we can write the partition 
function in the following form, 

Z = ZiD fvipViP^e-^'^'f'''f'^\ (4) 



where Zid is the partition function of the ID Hubbad model, and the effective action (free energy functional) S{4>, "0^) 
is given by 

oo 

S{^p,^^)^Si + J2Sn, (5) 

Si= J dxidx2 ^ ^l,,^ iXi)[t^ik^)S{xi - X2) - tl{k^)G<-^^ [Xi ~ X2M^k^ (X2), (6) 

Sn ^ lyn-l y^ / (^Xi...dX2nT±{k±i,...,k±2n)G'~"\xi,X2,.-.,X2n)i^l^k^i(^i)---^'y2^kj_2„{x2n), (7) 

2n 

T±{kj_i,...,k±2n) =Y\.^^^^^3^^J2^" k^j,0- (^) 



Here Xi = (a;||,i), G^^^{xi — X2) is the single-particle Green's function of the ID Hubbard model, and Gc is the 
connected 2n-body correlation functions of the ID Hubbard model. N± is the total number of the Hubbard chains. 
In this expression, all effects of electron correlations are incorporated into correlation functions of the ID Hubbard 
model. Then we can utilize the exact result of the ID Hubbard model. 

B. The case of infinitely long-range transverse hopping 



We, first, consider the limit of infinitely long-range transverse hopping, l± ^ 00. As was pointed out by Boies et 



'., in this limit, S'„>2 vanish. This is due to the fact that for l± —^ 00, tj_{k±) -^ {t±/\/D±) X^iJi ^k±j,o, and thus 



the factor 1/iV"""'^ = l/l]^~ ^ in eq.(R) remains even after the summation over k±. As a result, the action is given 
by 5*1, i.e. the non- interacting ferniion with the renormalized mass. The single-particle Green's function in this limit 
is given by. 

Although the ID Hubbard model is exactly-solvable, the exact expression for correlation functions is known only 
for the low-energy Luttingcr liquid statclljEil For simplicity, we_use the following approximated expression for the 
Green's function of the Luttinger liquid in the momentum spacejld 

"-^"'^ y/{v,k±e)iv,k±ey 

where Vc and Vs are the velocities of holon and spinon, respectively, and 6^ = \{Kc + ^ 2) with K^, the Luttinger 

liquid parameter. Eq is an ultra-violet cutoff energy for the Luttinger liquid state. Substituting this into eq.(P), one 
finds that away from half-filling, the quasiparticle pole of eq.(H) exists, and thus the system is the Fermi liquid. In 
this state, the Fermi surface is shifted by, 

Afc^ = f "':^tt^' )~sgnft^(fcj), (11) 



^VsVcEq" 

and the single-particle energy dispersion near the Fermi surface is given by. 



E,^2il-9,)^^^^^\k\=v\k\. (12) 



We have the finite quasiparticle weight. 



_ 2t^{k^F)E-^'^{vcAkp)^<>^ 
(Vc + Vs)AkF 



It is noted that the Fermi surface shift eq.(|llD is non-zero only at k±j = 0, and thus the change of the Fermi surface 

volume due to this shift vanishes for infinitely large D±. Namely, AVfs/Vfs = AkF/kp(2Tr)^^ > in the limit of 

D± — > 00. This is compatible with the Luttinger's theorem. 
In this special limit, the effective action is given by, 

-S*™/ = 51 / '^^\\'^^'^lkJ^ht)Zk^i9t - vd^^^Hak^{xi\,t). (14) 

k± •' 

This describes the Fermi liquid mean field solution (RPA solution) of the coupled chains. 

Using the these expressions, we can investigate correlation effects in the vicinity of MH transition. According to 
the Bethe ansatz exact solution, near the half-filling, Uc ~ (5 ((5 = 1 — n, n, electron density). Eq is roughly given by 
Vc{Tr — 2fcF)/2 ^ 5^ . It should be noted that the above expressions are appficable provided that the shift of the Fermi 
surface Akp is smaller than 7r(5/2. If this condition is not satisfied, we can not use the expression of the correlation 

function of the Luttinger liquid. Therefore our following argument is restricted to the case oi 5 > ctj^ , where c is 
a constant of order unity. However, for sufficiently small tj_, we can approach the vicinity of MH transition point. 



where strong correlation effects characterizing the metal-insulator transition manifest. Near the half- filling but for 
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5 > cij^ , the quasiparticle weight is given by, 



Thus the quasiparticle weight is much suppressed, and decreases toward zero, as the electron density approaches the 
half-filling. 

Now we investigate thermodynamic and transport properties in the vicinity of MH transition. Thermodynamic 
quantities are given by the sum of the contribution from the ID Hubbard model and the Fermi liquid state; i.e. 
Q = QiD + Qfl- In this paper, we concentrate on Qfl in order to see the correlation effects in the Fermi liquid 
state near the Mott transition. We can easily calculate the (5-dependence of the Drude weight Dq, the specific heat 
coefficient 70, and the compressibility kq of the Fermi liquid part in this non-interacting limit, 

Do = Y,vl5{^i-Ek)^5, (16) 

k 



70 



= Y,5{^Ji-Ek)^5-\ (17) 



KcO 



= ^z2^(M-i?fe)^<5-i^. (18) 



The behaviors of the Drude weight and the specific heat coefficient are similar to those of the ID Hubbard model. 
It is noted that £>o7o is independent of 5. This is easily understood since there are no vertex corrections in this 
non-interacting limit. 

For the ID Hubbard model, < 6'c < 1/8. Thus Kco shows a divergent behavior, as the electron filling approaches 
the half-filling. This is also due to the lack of vertex corrections. However, in contrast to the ID Hubbard model, 
where the mmpressibility behaves like k ~ 5^^, the exponent of Kco is non-universal depending on the strength of 
interaction.Q 

The single-particle density of states of electrons at the Fermi surface also shows a remarkable divergent behavior, 



n-2ee 



p{Q) = Y,^k5{ti-Ek)^5 ^(T^. (19) 

fc 

This result makes a shajp-jfoiitrast to the MH transition in the infinite-dimensional Hubbard model, where the density 
of state remains finitcJj't3L3 The divergent behavior of the density of states is also similar to the MH transition of 
the ID Hubbard model. However the exponent of eq.(|l^) is much modulated by electron correlation. 

III. FINITE-HOPPING RANGE CORRECTIONS 

In this section, we consider the l/^_L-corrections which generate 2n-body interaction as shown in Fig. 1. Then the 
model can not be solved exactly anymore. What we can do is to perform a perturbative expansion in terms of l/l±. 
First of all, we give a scaling argument and justify our perturbative approach in the next subsection. 

A. Renormalization group argument 

We consider correlation effects beyond the mean field solution given by Smf, eq. ([l4|) . We apply the scaling trans- 



formation which keeps Si = 5^/, eq.(14), invariant,E3 

X —+ Xx, t —* Ai, 

iJa{x,t) -. A"^,(Af,Af). (20) 

Then, the scaling dimension of ipaix,t) reads a — ^^^^ ■ Now we consider the scaling equation for the effective 
coupling of Sn {n > 2), gn ^ {t±/l±)'^"' . Applying the perturbative renormalization group argument, we obtain the 
scaling equation. 



^^"^ [3n-(n-l)i^i-7„].9„. (21) 



dlnA 

Here A is the dimensionless scaling parameter which goes to infinity in the low-energy limit, and 7ri > is the scaling 
dimension of Gc {xi^xof.^-, X2n)- Although we know the value of 7„ for the ID Hubbard model from the conformal 
field thoery argument ,E3l2I Gc ixi,X2, ...,X2„) in Sn>2 may show a deviation from the critical power law behavior in 
the long-distance limit owing to the deformation of the Fermi surface discussed above. In the conventional Ginzburg- 
Landau arguments of coupled chainsJHj it is often assumed that the coefficient of the fourth order term is a non-zero 
constant. This means 7„ = 0. However the concrete value of 7„ does not matter in the following argument. We 
just assume it is finite. Then we can easily see that for D± > (3n — 7„)/(n — 1), Sn with n > 2 is irrelevant in 
the low-energy state, and can be treated perturbatively. The existence of the critical dimension for the irrelevance 
of Sn>2 means nothing more than the stability of the RPA like mean field solution in the large dimensional system 
with the long-range inter-chain interaction, which is consistent with the common sense that a mean field solution is 
more stabilized in the presence of long-range interaction in spatial dimension larger than an upper critical dimension. 
In this paper, we consider the infinitely large D± which simplifies the pcrturbative calculation as well as ensures the 
controllable expansion. In this limit, the self-energy does not depend on the transverse momentum k±, f^^ the case 
of the infinitely large dimensional Hubbard model, where the self-energy is a completely local quantity.cJ'cS 

In the following, we investigate strong correlation effects by performing the pcrturbative calculation in l/l± Before 
doing that, we need to care about any possibility of spontaneous symmetry breaking such as magnetic transition 
which may be brought about by the interaction between quasiparticles. However, in this paper, we assume that no 
spontaneous symmetry breaking occurs before the MH transition, since we are concerned with correlation effects in 
the Fermi liquid state near the paramagnetic MH transition point. 

Before carrying out the pcrturbative calculation, we would like to make a comment about the Luttinger's theorem 
which holds for the case of finite ^j^. The change of the Fermi surface volume due to the Fermi surface shift, eq.(nl]), 
is given by AVps — 2 J2k ^kp. In the case of D± ^ oo, we have, 

/oo 
dt\t\^sgn{t)p_i_{t) = 0, (22) 

-OO 

where 



p^{t) ^ 5] -5(t - t^(fc^)) ^ -^e~^. (23) 

Thus the Fermi surface volume is preserved. This is consistent with the Luttinger's theorem. 

B. Perturbative calculation of the self-energy 

We evaluate the l//j_-corrections to the self-energy. We neglect the diagrams which just give a chemical poten- 
tial shift, as shown in Fig. 2(a),(b),(c). The diagram of order 0(5!) (Fig. 2(d)) is most important for the mass 
enhancement factor. The self-energy correction from this diagram is given by, 

^f{en)-gl5-^i^ne^ + ..., (24) 

The next order diagram which contributes to the mass enhancement factor is of order 0(g2)(Fig- 2(e)), since g^ <^ ^f 
as seen from eq.(pi|). It is evaluated as. 

Note that E^^^^ shows more singular (5-dependence than E'^'*^. Higher order corrections also give more singular 5- 
dependence. Thus we need to be careful with the applicability of our perturbative calculation. According to eq. (pl|) , 
g2 ~ 1/A^^"^. Then, the ratio of eq.(§|) to eq.(|4|) is sufficiently small provided that A > ((5-3^=/(i-^=))i/(-Di^. 
The applicability of the result of order 0(5!), (p4[), is restricted to this low-energy scale. In other words, as the 
electron density approaches the half-filling J — > 0, the energy range where eq.(p4) can be used becomes very narrower 
such as /S.E ^ i<;^((53^c/(i-ec)^i/(£'i-6) ^n the following, we restrict our argument to this energy range. 



C. Specific heat coefficient 

Using eq.(p3) we have the mass renormahzation factor, 



de 



l + SSc 

(5^crr^. (-26) 



The quasiparticle weight is much more suppressed compared to eq.(|ij). In contrast, the quasiparticle damping 
z^ImSfe(e) is still infinitesimally small in the low-energy state. Thus the Fermi liquid theoHf is applicable even in the 
vicinity of the MH transition. We obtain the ^-dependence of the specific heat coe0icient,c3 

7=yE'5(A'-^fc)-^'5-T-^. (27) 

k 

The specific heat coefficient shows the divergent behavior as 5 decreases. The mass enhancement is larger than that 
of the non-interacting limit obtained in §2. Note that the exponent of eq.(p7|) depends on the strength of interaction 
indicating the non-universal behavior. 

D. Single-particle density of states 

Up to the second order corrections in 92 , the single-particle density of states at the Fermi surface is given by, 

p(0)=^z,t<5(A*-i?,*)^5-w^. (28) 



The result is the same as that obtained in §2(eq.(19)). Higher-order corrections does not affect this result, as expected 
from the conventional Fermi liquid theory. 

E. Drude weight 

In the computation of the Drude weight, we take into account the effect of the velocity renormahzation, «^ = 
zl.dk{Ek -\- Rel]fe(0)) due to the self-energ« correction, eq.(^), as well as the vertex corrections to electric currents 
which is nothing but the back flow effect ,E3 

Jk = vl+Y,44''C{k,k')5{^i-El,)vl,. (29) 

k' 

Here r(fc, k') is the interaction between the quasiparticles with momenta k and k' . When we calculate the current 
vertex corrections, we generally need to pay attention to the charge conservation law. However the charge conservation 
is not trivial in our system. The non-perturbed state obtained in §2 is composed of the Fermi liquid part and the 
Luttinger liquid part. As a result, the weight of the Fermi liquid part may be transfered to the Luttinger liquid part 
by the interaction introduced in this section. The charge conservation law is not expressed in the closed form only 
in the Fermi liquid part. Thus, here, we just evaluate the vertex corrections up to the second order in g2, of which 
the diagrams are shown in Fig. 3. Note that each line in Fig. 3 is the renormalized single-particle Green's function 
with the quasiparticle weight z^, eq. (Pq). This renormahzation effect must be taken into account, because the use of 
the bare propagator eq. (0) for the calculation of the diagrams in Fig. 3 leads to unphysical enhancement of current 
vertex corrections for small 5. This treatment may be justified by the Ward-Takahashi identity, though we do not 
discuss this point furthermore because of the reason mentioned above. Then we have the Drude weight. 



D^e^Y.''lJk5{^l-El)^5'^^. (30) 



It is noted that the relation, D^ — (5-independent still holds approximately. This is because that the corrections due 
to the back flow effect is sub-leading in this systems in the vicinity of the MH transition and thus the leading behavior 
of the Drude weight is determined by the velocity of quasiparticles. 



F. Compressibility 

Here we will see that the effect of the interaction between quasiparticles affects the (5-dependence of the compress- 
ibility drastically. Using the dressed single-particle Green's function with the mass renormalization factor z^, as in 
the case of the Drude weight, we compute the lowest order diagram for the compressibility as shown in Fig. 4(a), 

Kc^6^. (31) 

Higher order diagrams shown in Fig. 4(b) give sub-dominant contributions for small S. 

In contrast to the result in §2, as (5 — > 0, Kc decreases toward zero, though we can not reach the half-filling point, 
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(5 = 0, because of the restriction ctj^'' < 6. In the conventional Fermi liquid theory, it is expected that vertex 
corrections due to electron-electron interaction suppress the compressibility. However, here, the lowest order bubble 
diagram Fig. 4(a) gives the strong suppression of the compressibility. This implies that the usual Fermi liquid effect of 
vertex corrections may be incorporated effectively into the mass renormalization factor z^ in our model. The drastic 
change of (5-dependence due to 1/Z^-corrections means that Sn>2 terms, eq.(|7|), are in a sense dangerously irrelevant 
for D^>Q. 

IV. EFFECTS OF DISORDER: A POSSIBILITY OF MOTT GLASS 

In this section, we study the effects of random potential in the vicinity of the Mott traosition. We are mainly 
concerned with a possible realization of the Mott glass state proposed by Orignac et alJlJ The Mott glass state 
is characterized by vanishing compressibility and the dynamical conductivity which obeys the Mott's law, namely, 
(j{oj) ~ w^ for small w. In this state, although a single-particle exciation has a charge gap, particle-hole excitation 
is gapless showing the characteristic behavior of Anderson localization. Thus the Mott glass state is regarded as an 
intermediate state between Mott insulator and Anderson localization. This state is similar ta Anderson localized state 
in the presence of long-range Coulomb interaction in the spatial dimension larger than 2.E2I However the Mott glass 
state is realized for the case of finite-range interaction and even for ID systems. To see a possible realization of this 
state, we add the disorder potential term to the Hamiltonian: 

i,3 

Here ^ij is a random potential. In the presence of disorder, we can derive the effective action in a simalr way to that 
applied for clean systems. Using the Hubbard-Stratonovich transformation, 

■^3,3' = ''J-ij' +"jj'?«ii 

^33"^ai,3^'yi,3' ~^ ^33' ^crij'^crij' ~^ ^33''^cri,j'^'^i,3' + ^ij' ^crij V'ffi j' j (33) 

we decouple the disorder term as well as the transverse hopping term. Then integrating out c, c' fields, we have, 

oo 

S{i;,i^^) = Si + Y,Sn, (34) 

n=2 

Si ^ dxldx2'^■^pl^.^{xl)[A{x^l, k±) 

-A(x^^-L,k^_)A{xp,k^)G^^\xi - X2Mak^{x2), (35) 

Sn = I Am-1 /. / dxi...dx2nB{xi\i,...,X\\2n;k±l,...,k±2n) 

xG'(")(a;i,a;2,...,X2„)V'iifc^i(a;i)---'0<T2„fei2„(a;2n), (36) 

2n 
B{xili,...,Xll2n;k±l,...,k±2n) = ]J A(X||j- , fc^j ) Jy-2„ ^^^ ^^j. (37) 

J ^ 

Since our model has the infinite D±, Anderson localization in the transverse direction does not occur. For this reason, 
we concentrate on effects of disorder only in the direction parallel to the chains, and consider randomness without 
spatial dependence in the transverse direction, namely. 



{C{x\\^)^{xl\k)) = D^S{x\l, - X||fc). (38) 

This assumption simplifies the following argument substantially. Moreover we consider sufficiently small disorder to 
stabilize the single-particle charge gap_]J_.is known that the Mott-Hubbard gap of the ID Hubbard model is not 
collapsed by sufficiently small disorder .tjlla In the following, we take into account only the backward scattering due 
to random potential and neglect the forward scattering which is not important for the realization of the Mott glass. 
Using the renormalization group argument similar to that developed in §3.1, we can easily show that Sn>2 of eg. Pq) 
are irrelevant in the low energy state. Thus we neglect the disorder effects of Sn>2, and take into account the lowest 
order corrections in 1//^ as done in §3. Then, we have the low energy effective action, 

Seff = / dx(hjil}'^^f,^^{x\\)[zl^^ {lj - v*pdx^^ ) - ^{x\\)]'ilJakj_Lj{x\\) + irrelevant interactions. (39) 

In the following, we neglect the irrelevant interactions of this action. As a result, the ipedel is equivalent to the ID 
disordered non-interacting fermion system, which has been studied extensively so far.E3'L2l However the important 
difference from the classical model is that the effects of electron correlation are iucprporated into the mass renoppal- 
ization factor z^. and the renormalized Fermi velocity Vp. Following Berezinskii,E2l and Abrikosov and RyzhkinJ22l we 
can caluculate the dynamical conductivity. The result is given by. 



aiu.)^:--f^[ln{^-^]] . (41) 



This expression shows the Mott's law, indicating the Anderson localized state. For small S, it behaves like, 

-Df- 

For the ID Hubbard model, < 9c < 1/8. Thus as (5 — > 0, it is expected that (7{uj) vanishes, though we can not reach 
6 = because of the reason mentioned in the previous sections. In this case, the Mott glass state may not realize at 
half-filling. However in the presence of a finite range interaction, such as next-nearest or next-next-nearest intearction, 
it is possible to be 9c = 1/4 for (5^-0. It is noted that even for 9c = 1/4, the single-particle Green's function, eq.(Q) 
still has a quasiparticle pole. Moreover the velocity of holon near the half-filling point behaves like Vc ^ S even in 
the presence of a finite-range interaction provided that the low-energy effective theory of the single chain is given by 
Sine-Gordon model. Thus the argument presented in the previous sections is applicable. Then, in this case, there is a 
possibility of the Mott glass state with (j{u!) ^ uj'^ for 6^0. This result is consistent with the physical intuition that 
the presence of a finite-range interaction reduces particle-hole excitation energy lower than a single-particle excitation 
energy, stabilizing the Mott glass state. Since our approach breaks down just at the half-fillig, the above result is not 
a direct evidence of the Mott glass state. However it implies a strong possibility of its realization at half-filling. 

V. SUMMARY AND DISCUSSION 

We have studied strong correlation effects of the Fermi liquid state near the MH transition in the quasi-one- 
dimensional Hubbard model. We have treated the intra-chain electron correlation exactly exploiting the exact results 
of the Bethe ansatz solution. The systematic perturbative expansion in terms of the inverse of the transverse hopping 
range has been used for the analysis of the inter-chain electron correlation effects. It has been shown that the 
precursor of the MH transition can be seen in the (5-dependence of the Drude weight, the specific heat coefficient, and 
the compressibility, and that the exponents which depend on the strength of electron-electron interaction characterize 
non-universal behaviors of the MH transition. We have also studied effets of disorder on this strongly correlated state. 
We otain the dynamical conductivity which obeys the Mott's law. The result implies a possible realization of the 
Mott glass state at half-filling. 

In real systems, the situation with a large l± hardly realizes. We would like discuss about to what extend our 
results are valid in systems with small l±. In principle, we can incorporate the finiteness of l± by the systematic 
expansion. As shown in §3, if we take into account the second-order self-energy corrections, we have the Fermi liquid 
state with a strongly renormalized quasi-particle weight, and then near half-filling, higher order corrections in l/l± 
are much suppressed, (see the calculation of the Drude weight and the compressibility in §3.) Thus it is expected 
that the physical picture of weakly-interacting quasi-particles with a strongly renormalized mass may be hold in the 



vicinity of the Mott transition even for the case with small l± , though the exponents of physical quantities should be 
quantatively changed from those obtained in §3. 
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As claimed in §2., the applicability of our results is restricted to S > ctj^ . Thus, unfortunately, we can not reach 
the half-filling point. We need the information for off-critical behaviors of correlation functions of the ID Hubbard 

2/3 

model to investigate the case of < (5 < cij_ . Recently, the substantial progress of the study of correlation functions 
of ID massive integrable models has been achieved.E3 It may be possible to apply the new technique for the case of 
< 5 < ct/ . We would like to address this issue in the near future. 
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FIG. 1. Some examples of multi-particle interaction. 
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FIG. 2. Diagrams for selfenergy corrections of order (a) 0{g2), (b) 0(53), (c) 0(54), (d) 0(02), (e) 0{g2)- 







FIG. 3. Diagrams for vertex corrections to the electric current up to the second order in 32. 
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FIG. 4. (a) Diagrams for compressibility without vertex corrections, (b) Some examples of diagrams for compressibility with 
vertex corrections 
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